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In 2012, the High Energy Stereoscopic System (H.E.S.S.) was expanded by a fifth telescope 
(CT5). With an effective mirror diameter of 28 m, CT5 is able to detect the Cherenkov light of 
very faint gamma-ray air showers, thereby significantly lowering the energy threshold of this tele¬ 
scope compared to the other four telescopes. Extracting as much information as possible from 
the recorded shower image is crucial for background rejection and to reach an energy threshold 
of a few tens of GeV. The camera of CT5 is conceived to register the time of the charge pulse 
maximum with respect to the beginning of the 16 ns integration window of each pixel. This infor¬ 
mation can be utilised to improve the event reconstruction. It also helps to reduce the background 
contamination at low energies. We present new techniques for background rejection based on 
CT5 timing information and evaluate their performance. 
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1. Introduction 

Within the last few decades, Imaging Atmospheric Cherenkov Telescopes (lACTs) have greatly 
enhanced the knowledge about the most violent processes in the universe. Using the Earth’s atmo¬ 
sphere as a calorimeter leads to a large detection area, which enables lACTs to study sources at 
y-ray energies from several tens of GeV up to around lOOTeV. The disadvantage of this approach is 
the dominant background which mainly corresponds to charged cosmic rays. Developing efficient 
mechanisms to suppress this component is thus crucial to improve the sensitivity of lACTs. 

Background suppression in ground-based gamma-ray astronomy is generally based on the 
characterisation of the shower image that is recorded by the camera of an lACT. However, several 
efforts were carried out to make use of the arrival time information of the Cherenkov photons as an 
additional, independent variable. In 2007, [1] exploited the timing information of the first MAGIC 
telescope using a template-based reconstruction method that resembles the one of [2]. Two years 
later the MAGIC collaboration presented several ways to improve the single-telescope performance 
by introducing several time-based variables [3]. 

H.E.S.S. is a system of five lACTs located in Namibia. Erom 2004 to 2012 (con'esponding to 
H.E.S.S. phase I), the system consisted of four identical lACTs (named CTl-4) with an effective 
miiTor diameter of 12m each, disposed on the corners of a square with 120m edge length. In 2012, 
the construction of an additional large-size (32 x 24m) telescope (CT5) in the centre of the array 
was completed, marking the beginning of H.E.S.S. phase II. 

The improved electronics of the CT5 camera can trigger events with a rate of up to 5 kHz. This 
high trigger rate combined with its greater resolution allows to have a performing monoscopic anal¬ 
ysis using only the CT5 images. In addition to the intensity, the 2048 camera pixels record temporal 
informations of the events, such as the time over a given threshold and the Time of Maximum. 

In this contribution, we present the status of the background rejection study using the Time 
of Maximum information of the CT5 camera. The variables commonly used by the Model recon¬ 
struction [2] are firstly reminded. Then the Time of Maximum, its calibration and its potential for 
background suppression are described. Einally, preliminary results of the approach as tested on 
Monte Carlo simulations and data of the AGN PKS 2155—304 are shown. 

2. Usual variables for background rejection 

In this section the principal methods of background rejection used within the H.E.S.S. mono¬ 
scopic Model reconstruction are summarized. 

2.1 Shape cuts 

After the shower reconstruction, some primary cuts, called shape cuts, are performed. They 
are not especially background rejection cuts but more similar to safety cuts allowing to reject non¬ 
physical or hardly analysable reconstructed shower images. They also contribute to considerably 
reduce the number of events to analyse. Here are presented the main shape cuts used by the Model 
analysis : 

• Minimum charge : the minimum number of photo-electrons accepted in the image. 
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• Nominal distance : distance between the centre of gravity of the shower and the camera 
center. This allow to only keep events for which the center of gravity is inside the camera. 


2.2 Goodness of fit 

We describe here the shower goodness which is the main hadron rejecter of the Model analysis 
and attests of the quality of the model adjustment. For more details on its derivation, see [2]. The 
goodness G of a.N pixel image is given by : 


G = 
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where Si and /r, are respectively the recorded signal in photoelectrons and the model derived from 
the reconstruction fit for a given pixel i. 

The likelihood L{s\lJ.) defined by fhe probabilify P{s\iJ.,ap,Gy) of having a signal s when 
knowing a model is given by : 


lnL(5|jU) = —2lnP{s\n,ap,ay 


= -21n 




n\J2n{ol + no^) 


■exp 


( {s-nf \ 

V 2(a| + na2)J 


( 2 . 2 ) 


where n is fhe number of phofoelecfrons, Op is fhe widfh fhe pedesfal and Oy is fhe widfh a fhe 
single phofoelecfron peak. 

Finally fhe average of fhe likelihood is : 


<lnL>|^ = j \nL{s\p) X P{s\p,Op,ay)As (2.3) 

2.3 Primary depth 

The primary depth is the atmospheric depth of the first Cherenkov emission of a shower. It is 
reconstructed during the model analysis and allows to discriminate gamma from hadrons but also 
electron showers. Indeed, the Cherenkov emission of electron showers is expected to begin higher 
in the atmosphere. 


2.4 Night Sky Background Goodness 

The Night Sky Background (NSB) Goodness is the most powerful cut in the case of mono- 
scopic analysis. Its main goal is to determine if the recorded image is compatible with an image 
which doesn’t contain Cherenkov photons, but mostly NSB photons which will enlarge the pixel 
distribution of the base line position of the electronic pedestal. 

The NSB Goodness (Gnsb) is very comparable in its construction to the goodness of fit, pre¬ 
sented in the section 2.2. The main difference is that we considerate a null model p: 

Y, (lnL(5|^ = 0)- < InL > |^=o) 

Gnsb = - -7=- (2.4) 
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Figure 1: Shower candidate in the CT5 camera. The intensity display, on the left, shows the shower position 
in the camera (high intensity pixels) and the distribution of pixels containing only NSB photons (intensity 
close to 0 photon-electron). On the right, the Time of Maximum display presents a clear difference between 
the NSB photons with a ToM randomly distributed and the Cherenkov photons from the shower presenting 
a ToM gradient within the camera. 

where the likelihood definitions are identical to the equations (2.2) and (2.3) by replacing /i by 0. 

In the monoscopic analysis, this variable rejects around 90% of the events allowing to produce 
well normalized map and high signal over background ratio and significance. Most of the rejected 
showers are low energy events, due to their low intensity which is more compatible with null model. 
Thus, cutting harder on the NSB goodness, directly increases the energy threshold of the analysis. 

In order to keep the low energy events without degrading the analysis, other alternatives should 
be considered. We present in the next sections one of them, using the timing information recorded 
by CT5. 

3. The temporal information in H.E.S.S. 

The 2048 pixels of the fifth H.E.S.S. camera record two temporal informations within the 16 
nanoseconds of the integration window. The first one is the Time over Threshold which represents 
the duration of the electronic pulse over a given charge threshold. The second one is the Time of 
Maximum which will be described into details in the following sections. 

3.1 Time of Maximum 

During the integration window, the temporal information of the pulse is sampled every nanosec¬ 
onds and each pixel records the temporal position of the pulse maximum, called Time of Maximum 
(ToM). For an atmospheric shower, the ToM is expected to be very close from pixel to pixel, when 
for the pixels containing only NSB, it will be randomly distributed. The difference between the 
Cherenkov and NSB photons ToM distribution is well illustrated in the figure 1 showing an exam¬ 
ple of a shower like event in the CT5 camera. 

3.2 Calibration of the Time of Maximum 

A calibration of the ToM over the camera is necessary because the 2048 pixels are distributed 
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over 256 redout cards, which are not perfectly temporally synchronised. It is effectuated using Flat 
Field runs. 

During these runs, a laser flash illuminates uniformly the entire CT5 camera. Due to the 
parabolic mount of the CT5 mirrors, the laser photons are expected to arrive in all pixels in the 
same time. This duration is well smaller than the electronics resolution. Comparing the Time of 
Maximum of the different pixels allows to reduce the disparities which can occur between them. 

After calibration the Root Mean Square (RMS) of the ToM distribution of all the pixels goes 
from 0.6 ns to 0.02 ns. 




Figure 2: Left; Distributions of for simulated protons (in blue) and gamma photons (in red) passing 
official H.E.S.S. cuts and Gnsb > 18. Right: Effective area for the three studied cut configurations at 18° 
zenith angle and source simulated at 1° off-axis from the camera center. 


4. Timing variables for background rejection 

The loosest official cut configuration applied in the Model analysis involves a NSB Goodness 
cut {Gnsb, see section 2.4) at a minimum of 24. Cutting at lower Gnsb, allows to gain a non neg¬ 
ligible amount of low energy events but also introduce large systematics effect on the background 
rejection, especially in the map productions. The timing variable presented in this section aims 
to lower the energy threshold of the H.E.S.S. experiment while keeping a controlled background 
subtraction. 

4.1 Temporal Night Sky Background goodness 

The conception of the temporal NSB goodness is very similar to its homonym in intensity 
(section 2.4). Therefore, the definition of the ToM NSB goodness (G^^^) is almost identical to the 
equation (2.4) : 
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Here we assume that the probability that a pixel has a Time of Maximum value ToM when 
a model T is expected is given by a Gaussian centered in T with a width a equal to the temporal 
resolution of a pixel obtained during the calibration. The likelihood is then given by : 

lnL(roM| t) = ln(27r) + In(a^) + ( 42 ) 

and the average likelihood by : 

< InL > It = ln(27r)+ ln(a^) + 1 (4.3) 

The goal is to compare the recorded image to an image that doesn’t contain Cherenkov photon. 
From a temporal point of view, the ToM of a null image is randomly distributed in the camera. The 
model T is thus picked randomly in a uniform distribution between 0 and 15 ns, in which is added 
a Gaussian temporal resolution. 

The loosest cut used in H.E.S.S. is obtained with the intensity NSB Goodness Gnsb > 24. In 
the next sections, we will lower it to Gnsb >18 and study the behaviour of the temporal NSB 
Goodness. This G^sb value has been taken to decrease sufficiently the energy threshold without 
increasing too much the map systematic uncertainties. 

4.2 Study of Monte Carlo simulation 

From the definition of the GJ^°^, its main purpose is to reject NSB triggered events and gamma 
candidate shower which are too much contaminated with NSB photons. In addition, this section 
will show that this variable also allows to improve the hadron suppression of the Model monoscopic 
reconstruction. The left of the figure 2 presenfs fhe GJJ’^ disfribufions for simulafed gamma and 
proton evenfs passing fhe usual H.E.S.S. cuf and Gnsb > 18. Protons presenf a higher rale lhan 
gamma al low Gjf^, meaning lhal fhe prolon shower liming dispersion is higher lhan gamma 
shower, as expected. Thus, a cuf al G]^^ > 1050, limil where fhe gamma rale becomes greater 
lhan Ihe proton rate, should improve Ihe background rejeclion, by suppressing more protons lhan 
gamma. 

The righl of Ihe figure 2 shows Ihe effeclive area graphics as funclion of Ihe simulated energy 
for fhe Ihree sludied configuralion : Ihe loosesl H.E.S.S cuf {Gnsb > 24), a drop in NSB Good¬ 
ness {Gnsb >18) and Ihis Iasi drop adding Ihe temporal NSB Goodness cuf found in Ihe previous 
paragraph {Gnsb > 18 > 1050). This plol shows lhal a rise of Ihe NSB Goodness direclly 

increases Ihe energy Ihreshold, while applying Ihe GJJ’^ cul has much less impacl on Ihe energy 
Ihreshold. 



Na 

S/B 

Ny 

Width of the OFF 
significance dislribulion 

Gnsb > 24 

18.4 

0.39 

1070 

1.0 

Gnsb >18 

17.2 

0.21 

1634 

1.3 

Gnsb > 18 + > 1050 

13.6 

0.24 

899 

1.1 


Table 1: Results of the PKS 2155-304 analysis for the three studied cut configurations. is the significance 
of the source, Ny is the number of excess events and S/B is the signal over background ratio. 
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Figure 3: Form left to right, the following cut configurations : Gnsb > 24, Gnsb >18 and G^sb > 
18 + Gjf^ > 1050. On top, significance distribution of the OFF events (in red) and excess events (in black). 
On bottom, excess maps of the PKS 2155-304 region. 

4.3 Application on PKS 2155-304 data 

The detection of the AGN PKS 2155-304 by H.E.S.S. was published in [4] based on the 2002 
data. Here the studied data set is constituted of 2014 data passing the standard H.E.S.S. run quality 
selection. It results in a total live time of 8.4 hours of CT5 data. 

The table 1 shows the results of the PKS 2155-304 analysis using the three sets of cuts de¬ 
termined in the precedent sections. Eowering the Gnsb cut from 24 to 18 increases the number of 
gamma candidates, but decreases significantly the signal over background ratio in the same time. 
The OPE event significance disfribufion, presenfed in fhe leff and middle of fhe figure 3, enlarges 
when decrease fhe Gnsb cuf value. The widfh of fhis disfribufion increases from a perfecf 1.0 a fo 

1.3 a. This behaviour has a direcf effecl on fhe excess map of fhe PKS 2155-304 region shown on 
leff and middle of fhe figure 3, which presenfs a nof expecfed bump on fhe fop leff of fhe source, 
when using fhe lowesf Gnsb cuf. 

Adding a G]^^ cuf af 1050 decreases significanlly fhe number of excess evenfs and fhe source 
significance, when fhe signal over background remain more or less sfable (fable 1). Buf fhis cuf 
improves fhe OPE evenf significance disfribufion and fhe map of fhe PKS 2155-304 region (righf 
of fhe figure 3). The widfh of fhe OPE significance disfribufion decreases from 1.3 a fo 1.1 a and 
fhe nof expecfed bump on fhe excess map disappears. 

Pinally, fhe figure 4 presenfs fhe reconsfrucfed energy disfribufion of fhe excess evenfs for fhe 
fhree sfudied cuf configurafions. The decreases of fhe energy fhreshold when lowering fhe Gnsb cuf 
from 24 fo 18 is once again verified. Adding fhe GJ^^ cuf lefs fhe energy disfribufion shape almost 
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unchanged, so the energy threshold doesn’t increase as expected from Monte Carlo simulation 
(see section 4.2). Compared to the Gnsb > 24 cut (official loosest cut of the Model analysis), the 
energy distribution using the ToM variable present more excess at low energy and less excess at 
high energy. This phenomenon can easily explain the decrease of the significance and number of 
excess events, since the cuts more high energy excess events than it gains at low energy. 

The time of arrival of the Cherenkov photons emitted by the highest energetic gamma is expected 
to be wider due to the larger bight of the electromagnetic shower, explaining the loss of the high 
energetic events when using a GJ^^ cut. 



Figure 4: Distribution of the reconstructed energy of the excess events 


5. Conclusion 

A new variable using the Time of Maximum recorded by the fifth H.E.S.S. camera has been 
introduced. This variable is constructed to reduced the contamination of shower events with NSB 
photons. Tests on Monte Carlo simulations and real data analysis of the PKS 2155-304 regions 
have shown that it can also be used for hadron rejection and allows to decrease significantly the 
energy threshold compared to the loosest official H.E.S.S. cuf configurafions, while keeping confrol 
on fhe background subfracfion. All resulfs presenfed in Ibis paper are preliminary. 
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